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The effect of gravitational radiation reaction on orbits around a spinning black hole is analyzed. 
Such orbits possess three constants of motion: i, e, and a, which correspond, in the Newtonian limit 
of the orbit being an ellipse, to the inclination angle of the orbital plane to the hole's equatorial 
plane, the eccentricity, and the semi-major axis length, respectively. First, it is argued that circular 
orbits (e = 0) remain circular under gravitational radiation reaction. Second, for elliptical orbits 
(removing the restriction of e = 0), the evolution of i, e, and a is computed to leading order in S 
(the magnitude of the spin angular momentum of the hole) and in M/a, where M is the mass of 
the black hole. As a decreases, l increases and e decreases. 

PACS numbers: 04.25.Nx, 04.30. Db 



I. INTRODUCTION 



The Earth-based Laser Interferometer Gravitational 
Wave Observatory-(LIGO-)VIRGO §,| network of 
gravitational wave detectors (which is now un- 
der construction) and the European Space Agency's 
planned space-based Laser Interferometer Space Antenna 
(LISA) pi will be used to search for and study the grav- 
itational waves from "particles" , such as neutron stars 
and small black holes, spiraling into massive black holes 
(mass M up to - 3OOM for LIGO/ VIRGO and up to 
~ 10 7 Mq for LISA). To search for the inspiral waves and 
extract the information they carry will require templates 
based on theoretical calculations of the emitted wave- 
forms, which in turn require a detailed understanding of 
how radiation reaction influences the orbital evolution. 

When the orbital plane of the particle is inclined to 
the equatorial plane of a spinning hole, only one method 
has been successfully implemented to deduce how ra- 
diation reaction influences the evolution of the orbit's 
"Carter constant" ||,[| , which governs the orbital shape 
and inclination angle. This method, which uses a "post- 
Newtonian" gravitational radiation reaction force, was 
described in a previous paper j(| , but there only applied 
to "circular orbits" (orbits of constant Boyer-Lindquist 
radial coordinate r) for simplicity. This follow-up paper 
has a two-fold purpose: First, in Sec. ||, we will argue 
that circular orbits remain circular under gravitational 
radiation reaction. Second, in Sec. |l|, we will compute 
the evolution of elliptical orbits under radiation reaction, 
but only to leading order in S, the magnitude of the spin 
angular momentum of the black hole, and leading order 
in M/a, where M is the black hole's mass and a is the 
size of the orbit, as defined more precisely below. (Here 
and throughout, units with G = c = 1 are used.) 
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II. EVOLUTION OF CIRCULAR ORBITS 



Several years ago, Ori H put forth the conjecture that 
circular orbits in the Kerr metric remain circular even 
under gravitational radiation reaction. Here, we will ar- 
gue in favor of the conjecture. We will start by review- 
ing some properties of elliptical and circular orbits in the 
Kerr metric. Then we will argue that a circular orbit and 
the reaction force acting on it have a type of reflection 
symmetry that ensures that the orbit remains circular 
under radiation reaction, in the limit of the particle's 
mass being small compared to the hole's mass. 

In the absence of gravitational radiation, the geodesic 
motion of a particle in orbit around a Kerr black hole 
is well-known and discussed, for example, in Sec. 33.6 of 
Rcf. The location of the particle can be described 
in Boyer-Lindquist coordinates r, 6, <f), and t. The orbit 
can be described by three constants of motion: the en- 
ergy E, the angular momentum along the hole's spin axis 
L z , and the Carter constant Q. The particle's rest mass 
/x can be counted as another constant of the motion. The 
energy E is defined as the relativistic energy of the par- 
ticle minus its rest mass, so that U E — /i" in the language 
of Ref. S corresponds to "E" here. We will restrict to 
bound orbits, that is E < and, as a consequence (see 
Ref. |), Q>0. 

An interesting feature of the Kerr metric in Boyer- 
Lindquist coordinates is the existence of nonequatorial, 
circular, geodesic orbits. These orbits are circular in 
the sense that the particle maintains a constant Boyer- 
Lindquist coordinate r; however, the plane of the cir- 
cular orbit is not fixed but rather precesses around the 
hole's spin axis. Such orbits exist and are stable for val- 
ues of E, L z , and Q that give R = 0, dR/dr = and 
d 2 R/dr 2 < 0, where R (see Eq. (33.33c) of Ref. g) is 
defined by 

R=[(E + At )(r 2 + S 2 /M 2 ) - L Z S/M] 2 
-(r 2 -2Mr + S 2 /M 2 ) 

x [ M V + (L z - Sn/M - SE/M) 2 + Q] . (1) 

For an arbitrary orbit with constants E, L Zl and Q, 
there might be some other energy E < E (E depends on 
L z and Q) such that, if the orbit had energy E(L Z ,Q) 
rather than E, the orbit would be circular and stable. In 
such a case, as an alternative set of constants to E, L z , 
and Q, the constants I, e, and a can be defined as follows: 



(2a) 



(Q + L 2 Z 



,1/2' 



l-e 2 ^|, (2b) 
Er 

a=— . (2c) 

Here f = r{L z , Q) is the radius of the circular orbit with 
constants E, L z , and Q. Note that a should not be con- 
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fused with the conventional notation for the spin of the 
black hole, which is S here. 

The positive root in Eq. ( pa| ) or in any other square 
root is always chosen. We choose the angle i to lie in 
the range < I < it, so that i < it/ 2 corresponds to 
an orbit co-rotating relative to the spin and l > tt/2 to 
counter-rotating. Also, e is chosen as nonnegative. 

This set of constants t, e, and a has the conceptual ad- 
vantage that in the Newtonian limit of large a, the orbit 
of the particle is an ellipse of eccentricity e and semi- 
major axis length a, on a plane with inclination angle 
i to the hole's equatorial plane. When not in the New- 
tonian limit, interpreting t, e, and a as the inclination 
angle, eccentricity, and semimajor axis length must be 
done with the caveat that since the orbit is not an el- 
lipse, then words such as "eccentricity" are subject to a 
modified interpretation and can be misleading. 

Even though the particle's motion is complicated when 
not in the Newtonian limit, some of the parameters that 
describe the particle's motion need not be specified. For 
example, we are not concerned with the value of <fi or 
t, because making a <f> or t translation does not change 
the physics in the axisymmetric, stationary Kerr metric. 
Another symmetry is that if the orbital motion is flipped 
over the hole's equatorial plane, i.e., 9{t) is replaced by 
ir — 9(t), the motion can be considered the same. All such 
</> and t translations and 9 reflections leave the shape of 
the orbit unchanged. 

We can think of the particle as undergoing oscillatory, 
coupled motion in the r and 9 directions. We define one 
orbital revolution to be one oscillation cycle as measured 
by the 9 motion. Given any chosen starting point of an 
orbital revolution with coordinate 9q, the revolution can 
be broken into two half revolutions, the first when the 
particle goes from 9q to ir — 9 half a 0-cycle later, and 
the second when the particle goes from it — 9 back to 9 
another half 6*-cycle later. [Because of the coupling of the 
r motion with the 9 motion, the 9 motion does not peak 
at the same extrema every cycle. Therefore, \tt/2—9q\ has 
to be chosen small enough that the orbit does indeed go 
through 7r — 9q and 9o in the following cycle. However, 
this is a very minor restriction for the rest of Sec. ||, 
where in proving that circular orbits stay circular, we 
only consider circular and almost circular orbits (we do 
not have to consider generally eccentric orbits since we 
know that a circular orbit cannot immediately become 
generally eccentric without first being slightly eccentric). 
In such case, the peaks of the 9 motion are almost the 
same every cycle.] 

Now we consider the effect of gravitational radiation 
reaction on an orbit. We assume that the rest mass \i 
is small enough for the adiabatic approximation to hold: 
the timescale of gravitational radiation reaction is much 
longer than any other timescale in the problem. Then 
the particle moves very nearly on a geodesic path char- 
acterized by the constants of motion i, e, and a; and only 
on a very long timescale (which varies like as /i — > 0, 
because the radiation reaction acceleration scales like fi) 
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is this motion substantially modified by gravitational ra- 
diation reaction. 

We now consider, for an orbit slowly inspiraling due to 
radiation reaction, an orbital revolution that satisfies the 
following condition, to which we give the name reflection 
symmetry: Consider the point on the orbit that is at the 
beginning of the orbital revolution. Denote by ro, 9o, ^o, 
#o, and 4>o the Boyer-Lindquist spatial coordinates of that 
point and their time derivatives. (Here and throughout, 
an overdot represents d/dt.) Then there are two other 
locations later on the path with coordinates 

r n = ro + irnr + h.o., (3a) 

9 n - tt/2 = (-1) 2 "(6» - tt/2 + /_m0) + h.o., (3b) 

r n = f + [inr + h.o., (3c) 

9 n = 6 + find + h.o., (3d) 

4> n = 4>o + ^n<t> + h-o., (3e) 

for rt = i (a half revolution after n = 0) and n = 1 (a full 
orbital revolution after n = 0). The functions with tildes 
are not functions of /i. The "h.o." terms are any terms 
that go to zero faster than \x as /i — * (higher order in \i 
than linear). 

Because of the initial conditions at the beginning of the 
first and second half revolutions (at n = and n = \, 
respectively), the shape of the first half revolution (the 
path connecting the n = and n = ^ locations) deviates 
from the shape of the second half revolution (the path 
connecting the n = | and n — 1 locations) by a path de- 
viation of order fi. Of course, these two paths also differ 
by a <f> translation, a t translation, and a reflection across 
the equatorial plane. But as we discussed above, these 
are unimportant differences because of the symmetries; 
the shapes of the paths are the same. 

Now that we have written Eqs. (||), we temporarily 
(for the remainder of this paragraph) go back to the case 
of no radiation reaction, i.e., we set to zero the y, terms 
and the h.o. terms in Eqs. (^). Clearly, a circular orbital 
revolution satisfies Eqs. (|j|) for any initial n — location 
chosen on the circular orbit. But could there be an ec- 
centric orbit which also satisfies Eqs. (||)? The answer 
is negative, as we shall now show. A slightly eccentric 
orbit (one with the value of e small enough that e 2 terms 
are negligible) would have the same 9(t) and <fi(t) mo- 
tion regardless of the value of e, but r — f would oscillate 
with an amplitude proportional to e. This can be verified 
from the Kerr-metric geodesic equations, Eqs. (33.32) of 
Ref. ||. In the Newtonian limit, the oscillation of r — f 
would be periodic with the same period that 9(t) has, 
but when not in the Newtonian limit the 9 and r oscilla- 
tions would have different periods. If an orbit were to be 
reflection symmetric, then r — f would have to have the 
same value when the orbital motion is at 9q as it does 
when it is at tt — 9q at the next value of n. This would 
require that either r — f oscillate at a frequency that is 
an even integer multiple of the 6 oscillation frequency, or 
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r — f have zero amplitude (a circular orbit). The former 
is never the case, as can be verified by numerically |§ ex- 
amining circular orbits in the Kerr metric over the space 
of possible physically acceptable values of S, L z and Q. 
The fact that r — f does not resonate with an even mul- 
tiple of the 9 frequency implies that a slightly eccentric 
orbit cannot be reflection symmetric. 

Now we shall return to the case of interest: that with 
gravitational radiation reaction. What precisely do we 
refer to when we discuss circular orbital revolutions, 
when the orbital revolution is not actually circular but 
is slowly inspiraling? A good, but not unique, defini- 
tion is one that agrees with the result in the case of no 
radiation reaction: We define that a circular orbital rev- 
olution is one that satisfies Eqs. (H), while an eccentric 
orbital revolution is one that does not (at least for slight 
eccentricity: as mentioned above, we are not considering 
generally eccentric orbits). An orbital revolution with 
weak radiation reaction is defined as circular if and only 
if it is reflection symmetric. 

We now consider starting with an initial orbital revo- 
lution that is circular, or equivalently, that is reflection 
symmetric, i.e., that satisfies conditions (|3|). For small /z, 
ignoring the h.o. (higher than /j.) corrections, we would 
expect that the third half revolution (the first half of the 
next orbital revolution) would have a shape that devi- 
ates from that of the second half revolution by the same 
amount as the shape of the second half revolution devi- 
ates from that of the first. We expect this, because from 
conditions (||) above, the initial conditions of the third 
half revolution differ from those of the second by the same 
amount (to linear order in fj,) as those of the second differ 
from those of the first; and the acceleration on the parti- 
cle should similarly be equally (also to linear order in fi) 
different between corresponding locations on the second 
and third half revolutions as between corresponding loca- 
tions on the first and second. The orbit remains circular 
for the additional half revolution. If there is any eccen- 
tricity added, it is in the h.o. terms, but in the /u — ► 
limit, this is ignorable compared to the shrinking of the 
orbit, which varies like fi (the terms involving tildes). 

We can repeat the above argument to get the shape 
of the fourth half revolution, as well as the fifth, sixth, 
etc. In fact, the argument can be repeated to any cho- 
sen number, n max , of orbital revolutions, as long as that 
chosen number does not go to infinity as fi — > 0; for if 
it did, then we would not be guaranteed that after the 
infinite number n max of orbits, the h.o. corrections of the 
above paragraph would be ignorable. For example, we 
could choose n max to be 100, but we could not choose it 
to be 100M//i. The orbit remains reflection symmetric 
(or equivalently, it remains circular) for n up to n max , 
where n increments by i. In other words, there is a lo- 
cation, with coordinates r n , 9 n r n , 9 n , and 4> n , satisfying 
Eqs. (|3|) for any n up to n max . 

The constants of motion E, L z , and Q (or equivalently, 
i, e, and a) evolve in such a way that in going from 
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n = to n = n max a circular orbit remains circular. By 
assigning new values of r*o, #o, ^o, 9q, and 4>o as the old 
r„ mMl #n max , r nm „, $n max , and 0n max , the argument can 
be repeated, over and over again. The rates of loss of 
E, L z , and Q will then continue at such a rate so as to 
maintain circularity. 

A more intuitive picture of why a circular orbit re- 
mains circular was provided by Ori [Q , who first pointed 
out that the incommensurability of the r and 9 periods 
is the key reason why the argument can be made with- 
out knowing the nature of the reaction force: Even if the 
radiation reaction were to take the bizarre form of some- 
body with a hammer hitting the particle every time the 
particle is at some value of 9, there would have to be 
another person across the equatorial plane at n — 9 with 
a hammer hitting the particle in a corresponding way, 
as dictated by the orbital symmetries. Since the r — f 
frequency is not an even multiple of the 9 frequency, the 
hammer hits cannot constructively interfere with each 
other and produce an eccentricity. 

If an orbit is circular, then just knowing the rates of 
change of E and L z (for example, by knowing the energy 
and angular momentum carried off in the gravitational 
waves) is enough to determine the full orbital evolution 
since the evolution of Q is constrained such that the con- 
ditions listed immediately before Eq. (|l|) are satisfied, for 
as long as the orbit itself is stable. 

III. LEADING ORDER EFFECT OF SPIN ON 
ECCENTRIC ORBITS 

We now wish to consider general, not just circular, 
orbits around a black hole. But in doing so, we restrict 
ourselves to only considering the leading order effect of 
spin. We will use the formalism of a radiation reaction 
force described in a previous paper Q and merely state 
how the method as described in that paper generalizes 
to orbits with eccentricity. 

When one is only interested in leading order in S and 
in M/r (or equivalently, M/a, in terms of orbit parame- 
ters), the effect of the hole's Kerr metric on the particle's 
motion can be substituted with a spin-orbit interaction 
in three dimensional flat-space. Let spherical polar coor- 
dinates r, 9, and 4>, centered on the black hole, be used 
to describe the location of the particle (these coordinates 
describe the relative separation of the two bodies), with 
the hole's spin along the polar axis. The Lagrangian 
(Ref. Eq. (4)) for the motion of the particle is given, 
to linear order in S but otherwise in solely Newtonian 
theory, by 

r , 2i& , 2 • 2/ fl \i2i , M M 2fiS sin 2 9 ■ 

L = —\r + r 9 +r sin (p)(b H (p. 

2 r r 

(4) 

To leading order in S and in M/r, the motion resulting 
from this Lagrangian is the same as in the Kerr metric. 
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The use of flat space coordinates, which ignores M/r cor- 
rections, is adequate to leading order. Using the same 
coordinate variable names r, 9, and (j> for these coordi- 
nates as for the Kerr metric's Boyer-Lindquist coordi- 
nates does not cause conflict and should not cause con- 
fusion. Alternatively, we can use Cartesian coordinates, 
x\ = rsin#cos(/>, X2 = r sin 9 sin (f>, and X3 = rcos9. 

The Lagrangian (|J) admits three constants of motion, 
called E, L z , and Q because they are the same constants 
as we have in the Kerr metric, to leading order in S and 
in M/r. The values of these constants are: 

E = £[f 2 + r 2 9 2 + r 2 sin 2 (0)</, 2 l - — , (5a) 
2 r 

t 2 • 2/m: 2[iSsin 2 9 

L z = /j,r sin {t))<p , (5b) 

Q + L 2 z =fi 2 r 4 [9 2 + sin 2 (9)^ 2 }~4fi 2 Sr sin 2 (9)<p. (5c) 

The combination Q + L 2 is a more natural constant to 
work with than Q. If S were equal to zero, then Q + L 2 
would be the square of the total angular momentum. 

The constants of motion 0, e, and a, when considered 
only to leading order in S and in M/r, are related to E, 
L Zl and Q by: 

cost = (6a) 

(Q + Ll) 1 ' 2 

1 - e 2 - _2 ^Q + £ «) (l + 4 ^ (6b) 

[j?M 2 \ + {Q + Ll) 2 )' (b j 

a ^wl 1 + 2 (QTWj- (6c) 

It is easy to verify these, by checking that the E and f 
that would make Eqs. (Q) give Eqs. (||) satisfy (at leading 
order in S and in M/a) the stable circular orbit con- 
straints listed immediately before Eq. ([j]). Note that 
Eqs. (^) are valid for arbitrary eccentricity e; they do 
not require e <g; 1. 

It is possible to express the instantaneous time deriva- 
tive of each constant of motion, dE/dt, dL z /dt, or 
d(Q + L 2 )/dt, as a function of r, xs, r, ±3, and the con- 
stants of motion; there is no occurrence of <f> (because of 
the axisymmetry) or <f> (as this is determined with L z , r, 
and 9 known) in any of the expressions. If S were zero 
then there could be no X3 dependence, rather only r de- 
pendence, since there is no physically preferred direction 
when spin is absent. Thus, an X3 or ±3 can only show up 
in a term that includes a factor of S. Because of this, to 
compute the time derivative of each constant of motion 
to Newtonian order plus the spin correction, xs(t) only 
needs to be known to Newtonian order, because the spin 
correction to x$ would be an S 2 term in the derivative 
of the constant of motion. On the other hand, the radial 
motion r(t) of the particle has to be known to Newtonian 
order plus the spin correction. The (f>(t) motion does not 
have to be known at all for computing the evolution of 
the constants of motion. 
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Let us, then, compute r and X3 to the necessary orders. 
One of the Euler-Lagrange equations yields 

The solution of this, in terms of a parameter ip, is 
(Q + Lp/{p?M) ( SL^M 



l + ecosV 1 V (Q + L l) 2 



(8) 



dt _(Q + Llf' 2 /{^M 2 ) f i + 6 SL £l lM_ . ( () ; 



# (1 + ecosT/.) 2 V (Q + ^z) 

In the Newtonian limit of S = 0, these are the equations 
for a Keplerian ellipse, with the true anomaly ip being 
the angle on the orbital plane of the particle relative to 
pcriastron as seen from the hole. 

To Newtonian order, X3 — r cos 8 can be expressed as 

Xs = r sint sin(^; + ipo). (10) 

Here, ipo is some constant that describes the orientation 
of the ellipse on the orbital plane. As seen from the hole, 
■00 is the angle between the direction of the periastron 
and the intersection of the equatorial and orbital planes. 
The orbital period, from periastron to periastron, is 

It happens that T, when written in this form, does not 
have an explicit S dependence. 

This motion we have just described is that in the ab- 
sence of gravitational radiation reaction; now we will 
compute the effect of the radiation reaction acceleration. 
We can take the equations for the rates of change of E, 
L z , and Q due to radiation reaction for a particle going 
around a more massive spinning body from Eqs. (10), 
(13), and (14) of Ref. 0. These equations give us formu- 
las for E, L z , and d(Q + L 2 z )/dt as functions of the dis- 
placement of the particle relative to the hole in Cartesian 
coordinates, Xk, and the relative velocity ik- There will 
also be higher order time derivatives of Xk (such as Xk , 
Xk, etc.), but these derivatives can be eliminated from 
the expressions for E, L z , and Q + L\ with the aid of 
the Euler-Lagrange equations [derived from (|j) when ex- 
pressed in Cartesian coordinates — note that repeated in- 
dices are summed over 1,2,3]: 

M ( 4 f L z 

(12) 

The time evolution of each constant of motion can 
thereby be expressed in terms of r, r, X3, ±3, and the con- 
stants of motion. The trajectory (p|)-(pd|) can be inserted 
into these expressions, and then time averaged using 



8 



E) = - 
T 



2tt 



, , dt ■ 



(13) 



and similarly for L z and Q + L 2 . The result is 
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to the other set of constants: 
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Equations (|14|) agree (after trivial conversions of nota- 
tion) with previous results: Eqs. (15) of Ref. || and (the 
first line of) Eq. (3.14) of Ref. pj, each of which is a 
special case of Eqs. (|l4|). 

In most cases, the terms with the cos(2?/>o) can be 
dropped because they average to zero; to see when this 
can be done, consider the following: The Newtonian ap- 
proximation to the motion is that the particle travels in 
an ellipse. The first correction to this motion is, as Ein- 
stein computed for Mercury, that the periastron position 
of the ellipse shifts on a timescale of 



T, 



prcc 



M(M/a)- 5/2 (l-e 2 ). 



(16) 



The radiation reaction timescale for terms that involve 
ifto , as computed by evaluating (sin t)j (L}^ terms, is 



Trad - M 



-11/2 



M 2 



(1-e 2 ) 



2 \ 1 



(17) 



(There are also factors of order unity that involve e which 
were ignored. If T ra( j were computed differently, for ex- 
ample by evaluating L Z /{L Z )^, terms, it would contain 
factors of i as well.) 

In the Newtonian limit, is fixed, but with the pe- 
riastron precession, -0 O changes slightly after each orbit, 



by a post-Newtonian correction that was ignorable until 
now: When T ra( j > T picc , the cos(2^> ) in Eqs. (14b) and 
( |l5a| ) averages to zero, and the terms with that factor can 
be dropped. For extremely eccentric orbits, T ra d might 
not be much greater than T prec , so the ipo terms must 
be kept. In all other respects, the periastron precession 
can be ignored because it just gives terms higher order 
in M/a (terms we have neglected). The only effect of 
the precession, to which our analysis is sensitive, is the 
averaging a way o f ipg in the case that T ra( j 3> T proc . 

From Eq. (15a), it is clear that the angle i changes such 
as to become antialigned with the spin. In Ref. ||, this 
conclusion was reached for circular orbits; finite eccen- 
tricity does not change, but only enhances, this result. 
However, the statement that "the inclination angle an- 
tialigns with the spin" is subject to the warning that we 
mentioned above when introducing l: With the orbit not 
confined to a fixed plane, the angle i is not the only way 
we could de fine "inclination angle" ||. 

Equation (15c) has two important consequences: First, 
to leading order, orbits tend to circularize, as is a well- 
known fact. Second, if an orbit is circular, then e = and 
(e) = 0, so the orbit remains circular. This is expected, 
since this is the leading order limit of the general result 
in Sec. || 

The above analysis is just one step in a general pro- 
gram for understanding the effects of radiation reaction 
on orbiting, spinning bodies. Future steps in this pro- 
gram include: generalizing the analysis to an arbitrary 
mass ratio n/M and to the case of both masses having 
spin, extending the analysis to higher order in M/r and 
in S, and achieving a similar calculation of the orbital 
evolution in the fully relativistic Kerr metric. 
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